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Abstract 

A general phenomenological theory is presented for the phase behavior of ferromag- 
netic superconductors with spin-triplet electron Cooper pairing. The theory accounts 
in detail for the temperature-pressure phase diagram of ZrZn2, while the main features 
of the diagram for UGe2 are also described. Quantitative criteria are deduced for the 
U-type (type I) and Zr-type (type II) unconventional ferromagnetic superconductors 
with spin-triplet Cooper electron pairing. Some basic properties of quantum phase 
transitions are also elucidated. 



The remarkable coexistence of itinerant ferromagnetism and unconventional (spin-triplet) 
superconductivity at low temperatures (T < 1 K) was discovered experimentally in the 
intermetallic compounds UGe 2 [H [21 E], ZrZn 2 [I], and URhGe Other metallic com- 
pounds, such as UCoGe 0, [7] and Ulr [HE], were also found to be spin-triplet ferromagnetic 
superconductors. In ZrZn 2 , URhGe, and UCoGe, the mixed phase of coexistence of ferro- 
magnetism and unconventional superconductivity (labeled the FS phase) occurs over a wide 
range of pressure (i.e., from ambient pressure P ~ 1 bar up to 10 kbar). By contrast, in 
other compounds (e.g., UGe 2 and Ulr) this FS phase is found only in the high-pressure part 
(P ~ 10 kbar) of the T — P phase diagram. 

Another feature of the above compounds is that the FS phase occurs only in the ferro- 
magnetic phase domain of the T — P diagram. Specifically, at equilibrium and a given P, 
the temperature Tp{P) of the normal-to-ferromagnetic phase (or N-FM) transition is never 
lower than the temperature Tps{P) of the ferromagnetic-to-FS phase (or FM-FS) transition. 
This is consistent with the point of view that the superconductivity in these compounds is 
triggered by the spontaneous magnetization M, by analogy with the well-known triggering 
of the superfluid phase Ai in 3 He at mK temperatures by the external magnetic field H . 
This helium-analogy has been used in some theoretical studies (see, e.g., pHl dU [121 02]), 
where Ginzburg-Landau (GL) free energy terms to describe the FS phase were derived by 
symmetry arguments. 
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For the spin-triplet ferromagnetic superconductors the trigger mechanism was recently exam- 
ined in detail [13]. The main system properties are affected by a term in the GL expansion of 
the form M\tl>\ 2 , which represents the interaction of M = {Mj;j = 1, 2, 3} with the complex 
superconducting vector field if) = {ipj}. Specifically, this term triggers -0 7^ for certain 
T and P values. An analogous trigger mechanism is familiar in the context of improper 
ferroelectrics |14j . 

A crucial consideration in this work is the nonzero magnetic moment of the spin-triplet 
Cooper pairs of the electrons. While the spin-singlet Cooper pairs have net spin zero and 
are quite sensitive to the magnitude of the magnetic induction B, the spin-triplet pairs are 
known to be robust with respect to relatively large B. The phenomena of spin-triplet super- 
conductivity and itinerant ferromagnetism are both due to the same electron bands of the 
compounds: the /-band electrons in uranium-based compounds and the (i-band electrons in 
ZrZn2- However, the microscopic band theory of magnetism and superconductivity in non- 
Fermi liquids of strongly interacting heavy electrons is either too complex or insufficiently 
developed to describe the complicated behavior in itinerant ferromagnetic compounds. Con- 
sequently, several authors (see [TQl dH [321 ESI EE]) have explored a phenomenological de- 
scription within self-consistent mean field theory, and we build on a similar approach here. 

In this Letter, by focusing on ZrZn 2 and UGe2 with their contrasting types of behavior, 
we show that the T — P phase diagrams of spin-triplet ferromagnetic superconductors can 
be successfully described starting from the general GL free energy F(ift, M) established in 
Refs. [TUl [TTT [T2l [13]. The present phenomenological approach includes both mean-field 
and spin-fluctuation theory (SFT), as in [12], considerations. We propose a simple, yet 
comprehensive, modeling of the P dependence of the free energy parameters, from which it 
is shown that the phase diagram of ZrZn2 is obtained in good quantitative agreement with 
the experimental data [4j. Further, the main features [1] of the T — P diagram of UGe2 are 
also well-described within our approach. The theory is capable of outlining several different 
possible topologies for the T — P phase diagram, depending on the GL parameters of the 
material which can be chosen in accordance with experiment. Quantitative criteria emerge 
for two distinct types of behavior for unconventional ferromagnetic superconductors, which 
we label Zr-type and U-type. Further possible applications are to URhGe, UCoGe, and 
Ulr. Our results address questions regarding the order of the quantum phase transitions at 
ultra-low and zero temperatures. They also pose intriguing questions pointing to further 
experimental investigations of (e.g.) the detailed structure of the phase diagrams in the 
high-P/low-T region. 

Following Ref. [13] the free energy per unit volume, F/V = f(tj>, M), can be written in the 
form 



m M) = a^| 2 +||V| 4 +y |V> 2 i 2 +! E I^| 4 +«/M 2 +|m 4 +z 70 M.(^x^)+5M 2 |V| 2 - 

(1) 

The material parameters satisfy, as in [13], b s > 0, bf > 0, a s = a s (T — T s ), and = 
azf[T n — Tj? (P)], where n = 1 gives the standard form of a/, and n = 2 applies for SFT [T6] . 
The terms proportional to u s and v s describe, respectively, the anisotropy of the spin-triplet 
electron Cooper pairs and the crystal anisotropy. Next, 70 ~ J (with J > the ferromagnetic 
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exchange constant) and 5 > are parameters of the i^-M interaction terms. Previous mean- 
field studies have shown that the anisotropies represented by the u s and v s terms in Eq. ([I]) 
slightly perturb the size and shape of the stability domains of the phases, while similar effects 
can be achieved by varying the b s factor in the b s \il)\ A term. For these reasons, in the present 
analysis we ignore the anisotropy terms, setting u s = v s = 0, and consider b s = b > as an 
effective parameter. Then, without loss of generality, we are free to choose the magnetization 
vector to have the form M = (0, 0, M). 

A convenient dimensionless free energy can now be defined by / = //(6/Mq), where M = 
[afTfy/bf] 1 / 2 > is the value of M corresponding to the pure magnetic subsystem {rjj = 0) 
at T = P = and Ty = Tf(0). On scaling the order parameters as m = M/M and 
cp = ^/[(bf/by/^Mo] we obtain 

~ d) Tft 

f = r<p 2 + y + tm 2 + — + 27m0i0 2 sin# + 7l m 2 2 , (2) 

where <pj = \<fj\, <p = \(f\, and 6 is the phase angle between the complex <p 2 and The 
dimensionless constants are t = [T n — Tj(P)}, r = k(T — T s ) with k = a s b x j 2 j a^b x l 2 T^ x , 
7 = 7o/[«/^/o^] 1//2 ; an d 7i = bj^bbf) 1 ! 2 . The reduced temperatures are T = T/T/o, Tf(P) = 
T f (P)/T f0 ,f s (P)=T s (P)/T f0 . 

The analysis involves making simple assumptions for the P dependence of the t, r, 7, and 71 
parameters in Eq. fl2]). Specifically, we assume that only Tf has a significant P dependence, 
described by f f (P) = (1 - P) 1 ^, where P = P/P and P is a characteristic pressure 
deduced later. In ZrZn 2 and UGe2 the Pq values are very close to the critical pressure P c at 
which both the ferromagnetic and superconducting orders vanish, but in other systems this 
is not necessarily the case. As we will discuss, the nonlinearity (n = 2) of Tf(P) in ZrZn 2 
and UGe2 is relevant at relatively high P, at which the N-FM transition temperature Tp(P) 
may not coincide with Tf(P). 

The simplified model in Eq. (J2J) is capable of describing the main thermodynamic properties 
of spin-triplet ferromagnetic superconductors. There are three stable phases: (i) the normal 
(N) phase, given by <fi = m = 0; (ii) the pure ferromagnetic (FM) phase, given by m = 
(—t) l l 2 > 0, 4> = 0; and (iii) the FS phase, given by 4>\ = 4>\ = (7m — r — 'j 1 m 2 )/2, <p 3 = 0, 
where sin^ = — 1 and m satisfies 



3 9 / 7 2 \ 7r 

'1 — 7 2 )m 3 + -77xm 2 + I t — - — 7ir j m + — = 0. (3) 



2 ' ' V 2 / 2 

We note that FS is a two-domain phase as discussed in Ref. [121 EE]. Although Eq. (j3J) 
is complicated, some analytical results follow, e.g., we find that the second order phase 
transition line Tps(-P) separating the FM and FS phases is the solution of 

f FS {P) = f s + —t(Tps) + 2 H(T F5 )] 1/2 - (4) 

K K 

Under certain conditions, the Tps{P) curve has a maximum at T m = T s + (7 2 /4k7!) with 
pressure P m found by solving t(T m ,P m ) = — (7 2 /47 2 ). Examples will be given later, but 
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generally this curve extends from ambient P up to a tri-critical point labeled B, with coor- 
dinates (Pb,Tb), where the FM-FS phase transition occurs at a straight line of first order 
transition up to a critical-end point C. The lines of all three phase transitions (N-FM, N-FS, 
and FM-FS) terminate at C. For P > Pc the FM-FS phase transition occurs on a rather 
flat, smooth line of equilibrium transition of first order up to a second trier itical point A 
with Pa ~ Po and Ta ~ 0. Finally, the third transition line terminating at C describes 
the second order phase transition N-FM. The temperatures at the three multi-critical points 
correspond to f A = f„ f B = f s + 7 2 (2 + 7i )/4k(1 + 7i ) 2 , and f c = f s + 7 2 /4/<l + 7l ), 
while the P values can be deduced from the previous equations. These results are valid 
whenever P/(P) > T S (P), which excludes any pure superconducting phase (ip 7^ 0,m = 0) 
in accord with the available experimental data. Note that, for any set of material parameters, 
T A <T C <T B < T m and P m < P B < P c . 

A calculation of the T — P diagram from Eq. (|2j) for any material requires some knowl- 
edge of Po, Tfo, T s , k 7, and 7 i. The temperature T/o can be obtained directly from the 
experimental phase diagrams. The model pressure Po is either identical to or very close to 
the critical pressure P c at which the N-FM phase transition line terminates at T ~ 0. The 
characteristic temperature T s of the generic superconducting transition is not available from 
the experiments and thus has to be estimated using general consistency arguments. For 
T f (P) > T S (P) we must have T S (P) = at P > P c , where T f (P) < 0. For < P < P , 
T s < Tc and therefore for cases where Tc is too small to be observed experimentally, T s 
can be ignored. For systems where Tc is measurable this argument does not hold. This is 
likely to happen for T s > (for T s < 0, T c is very small). However, in such cases, pure 
superconducting phase should be observable. To date there are no experimental results re- 
ported for such a feature in ZrZn2 or UGe2, and thus we can put T s = 0. We remark that 
negative values of T s are possible, and they describe a phase diagram topology in which the 
FM-FS transition line terminates at T = for P < P c . This might be of relevance for other 
compounds, e.g., URhGe. 

Typically, additional features of the experimental phase diagram must be utilized. For 
example, in ZrZn 2 these are the observed values of Tps(0) and the slope po = [dTFs(P)/dP]o 
at P = 0. For UGe2 one may use T m , P m , and Po c , where the last quantity denotes the other 
solution (below P c ) of Tps{P) = 0. The ratios 7 /k and 7 i/k can be deduced using Eq. (jlj) 
and the expressions for T m , P m , and p , while k is chosen by requiring a suitable value of 

Experiments [4J for ZrZn 2 indicate T f0 = 28.5 K, T FS (0) = 0.29 K, P ~ P c = 21 kbar, and 
Tp(P) ~ Tf(P) is almost a straight line, so n = 1 describes the P-dependence. The slope for 
Tps(P) at P = is harder to estimate; its magnitude should not exceed Tf /P c w 0.014 on 
the basis of a straight-line assumption, implying —0.014 < p < 0. However, this ignores the 
effect of a maximum, although it is unclear experimentally in ZrZn 2 , at (T m ,P m ). If such 
a maximum were at P = we would have p — 0, whereas a maximum with T m ~ Tps(0) 
and P m Po provides us with an estimated range < po < 0.005. The choice po = gives 
7//C ~ 0.02 and 7 i/k ~ 0.01, but similar values hold for any |po| < 0.003. The multi-critical 
points A and C cannot be distinguished experimentally. Since the experimental accuracy [1] 
is less than ~ 25 mK in the high-P domain (P ~ 20 — 21 kbar), we suppose that Tc ~ 10 
mK, which corresponds to k ~ 10. We employed these parameters to calculate the T — P 
diagram using p = and 0.003. The differences obtained in these two cases are negligible, 
with both phase diagrams being in excellent agreement with experiment. 
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Figure 1: T — P diagram of ZrZn 2 calculated taking T s = 0, p = 0.003 K/kbar, T/ = 28.5 
K, Pq = 21 kbar, k — 10, 7//? = 2ji/k ~ 0.2. The low-T domain of the FS phase is seen 
more clearly in the following figure. 



The latter value is used in Fig. 1, which gives Pa ~ -P c — 21.10 kbar, Pg = 20.68 kbar, 
P c = 20.99 kbar, T A = T F (P C ) = T FS (P C ) = K, T B = 0.077 K, T c = 0.038 K, and 
Tfs(Q) — 0.285 K. The low-T region is seen in more detail in Fig. 2, where the A, B, C 
points are shown and the order of the FM-FS phase transition changes from second to first 
order around the critical end-point C. The Tp${P) curve has a maximum at P m = 6.915 
kbar and T m = 0.301 K. These results account well for the main features of the experimental 
behavior [4], including the claimed change in the order of the FM-FS phase transition at 
relatively high P. Within the present model the N-FM transition is of second order up to 
Pc ~ P c - Moreover, if the experiments are reliable in their indication of a first order N-FM 
transition at much lower P values, the theory can accommodate this by a change of sign of 
bf, leading to a new tri critical point located at a distinct P tr < Pc on the N-FM transition 
line. Since Tc > a direct N-FS phase transition of first order is predicted in accord with 
conclusions from de Haas-van Alphen experiments [T7j and some theoretical studies [18]. 
Such a transition may not occur in other cases where Tc = 0. In SFT (n = 2) the diagram 
topology remains the same but points B and C are slightly shifted to higher P (typically by 
about 1 bar). 

The experimental data [H [21 [3] for UGe2 indicate T/o = 52 K, P c — 1.6 GPa (= 16 kbar), 
T m = 0.75 K, P m « 1.15 GPa, and Pq c ~ 1-05 GPa. Using again the variant n = 1 for 
Tf(P) and the above values for T m and Pq c we obtain j/k ~ 0.098 and Ji/k ~ 0.168. The 
temperature Tc ~ 0.1 K corresponds to k ~ 5. Using these, together with T s = 0, leads to 
the T — P diagram in Fig. 3, showing only the low-T region of interest. We obtain = K, 
T B = 0.481 K, T c = 0.301 K, P A = 1.72 GPa, P B = 1.56 GPa, and P c = 1.59 GPa. There is 
agreement with the main experimental findings, although P m corresponding to the maximum 
(found at ~ 1.44 GPa in Fig. 3) is about 0.3 GPa higher than suggested experimentally. If 
the experimental plots are accurate in this respect, this difference may be attributable to 
the so-called (T x ) meta-magnetic phase transition in UGe2, which is related to an abrupt 
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Figure 2: Detail of Fig. 1 with expanded temperature scale. 




Figure 3: Low-T part of the T — P diagram of UGe2 calculated taking T s = 0, Ty = 52 K 
P = 1.6 GPa, T m = 0.75 K, P Qc = 1.05 GPa (k = 5, 7//C = 0.098, and 71/K = 0.168). 



change of the magnetization in the vicinity of (T m ,P m ). Thus, one may suppose that the 
meta-magnetic effects, which are outside the scope of our current model, significantly affect 
the shape of the Tfs(P) curve by lowering P m (along with Pb and Pc)- It is possible to 
achieve a lower P m value (while leaving T m unchanged) , but this has the undesirable effect of 
modifying P c0 to a value that disagrees with experiment. In SFT (n = 2) the multi-critical 
points are located at slightly higher P (by about 0.01 GPa), as for ZrZn 2 . 

The estimates for UGe2 imply jiK ~ 1.9, so the condition for Tps(P) to have a maximum 
found from Eq. is satisfied. As we discussed for ZrZn 2 , the location of this maximum 
can be hard to fix accurately in experiments. However, P c q can be more easily distinguished, 
as in the UGe2 case. Then we have a well-established quantum (zero-temperature) phase 
transition of second order, i.e., a quantum critical point [19]. From Eq. (j3J) the existence 
of this type of solution in systems with T s = (as UGe 2 ) occurs for 7 < 71. Such systems 
(which we label as U-type) are essentially different from those such as ZrZn 2 where 71 < 7 
and hence Tps(0) > 0. In this latter case (Zr-type compounds) a maximum T m > may 
sometimes occur, as discussed earlier. We note that the ratio 7/71 reflects a balance effect 
between the two tp-M interactions. When the trigger interaction (typified by 7) prevails, 
the Zr-type behavior is found where superconductivity exists at P = 0. The same ratio can 
be expressed as 7o/5M , which emphasizes that the ground state value of the magnetization 
at P = is also relevant. In general, depending on the ratio of the interaction parameters 7 
and 71, the ferromagnetic superconductors with spin-triplet Cooper fermion pairing can be 
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of two types: type I (7 < 71) and type II (7 > 71). The two types are distinguished in their 
thermodynamic properties. 

The quantum phase transition near P c is of first order. Depending on the system properties, 
Tc can be either positive (when a direct N-FS first order transition is possible), zero, or 
negative (when the FM-FS and N-FM phase transition lines terminate at different zero- 
temperature phase transition points). The last two cases correspond to T s < 0. All these 
cases are possible in Zr- and U-type compounds. The zero temperature transition at P c o is 
found to be a quantum critical point, whereas the zero-temperature phase transition at P c is 
of first order. As noted, the latter splits into two first order phase transitions. This classical 
picture may be changed through quantum fluctuations [12J. An investigation [2QJ EE] of 
the quantum fluctuations and the quantum dimensional crossover by renormalization group 
methods revealed a fluctuation change in the order of this transition to a continuous phase 
transition belonging to an entirely new class of universality. However, this option exists 
only for magnetically isotropic order (Heisenberg symmetry) and is unlikely to apply in the 
known spin-triplet ferromagnetic superconductors. 

Even in its simplified form, this theory has been shown to be capable of accounting for a 
wide variety of experimental behavior. A natural extension to the theory is to add a M & 
term which provides a formalism to investigate possible metamagnetic phase transitions [22] 
and extend some first order phase transition lines. Another modification of this theory, with 
regard to applications to other compounds, is to include a P dependence for some of the 
other GL parameters. 

Finally, let us emphasize that our results are based on a general thermodynamic analysis 
in which the surface and bulk phases are treated on the same footing. For this reason, our 
results do not contradict to experiments [23J showing a lack of bulk superconductivity in 
ZrZn2 but the occurrence of a surface FS phase at surfaces with higher Zr content than that 
in ZrZn2. 
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